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We consider a directed polymer of length L in a random medium of space dimension d = 1, 2, 3. 
The statistics of low energy excitations as a function of their size I is numerically evaluated. 
These excitations can be divided into bulk and boundary excitations, with respective densities 
p b L ulk {E = 0,1) and p b ° mdar V(E = ,Z). We find that both densities follow the scaling behavior 
pb uik, boundary^ _ q ^ _ L - !-««, R bmk, boundary ^ x = y L ^ where ^ is th<J exponent governing the 

energy fluctuations at zero temperature (with the well-known exact value 9i — 1/3 in one dimen- 
sion).' In the limit x = l/L -> 0, both scaling functions R bulk (x) and R b °^ d ^y( x ) behave as 
R buik,boundary^ _ x -i-e d ^ i ea ding to the droplet power law p *»tk, boundary^ = ,;) ~ r 1 "^ in 
the regime 1 C l< L. Beyond their common singularity near x — > 0, the two scaling functions 
R buik, boundary^ are very dif f erent . w h e reas R bulk (x) decays monotonically for < x < 1, the 
function p l boundal v (x) first decays for < a; < x m in, then grows for x m i n < x < 1, and finally 
presents a power law singularity R boundar v ( x ) ^ _ a;) _ °"<i near x — > 1. The density of excitations 
of length I — L accordingly decays as pj°"™ ary (E = 0, 1 = L) ~ L _Ad where Ad = 1 + #d — <7d. We 
obtain Ai ~ 0.67, A2 — 0.53 and A3 ~ 0.39, suggesting the possible relation Ad = 29d- 

I. INTRODUCTION 

The model of a directed polymer in a random medium has attracted a lot of attention in the last twenty years for 
two main reasons. On the one hand, it is directly related to non-equilibrium growth models [|; on the other hand, it 
plays the role of a 'baby spin glass' model in the field of disordered systems [H 13, I Sfl J M - This model presents a low 
temperature disorder dominated phase, where the order parameter is an 'overlap' |2|, Ui |6l h| • This low temperature 
phase displays an extreme sensitivity with respect to temperature or disorder changes |ol ISl Wi ITfl Hdl Hp . fl^ | . and 
aging properties for the dynamics [14| . In finite dimensions, a scaling droplet theory was proposed |a, [l5j , in direct 
correspondence with the droplet theory of spin glasses |l6j |. whereas in the mean-field version of the model on the 
Cayley, a freezing transition very similar to the one occurring in the Random Energy Model was found |2J . The phase 
diagram as a function of space dimension d is the following [IJ : for d > 2, there exists a phase transition between 
the low temperature disorder-dominated phase and a free phase at high temperature (T?l ITq . This phase transition 
has been studied numerically in d — 3 [HI. l2fjj , exactly on a Cayley tree and on hierarchical lattice [2lJ . On 
the contrary, in dimension d < 2, there is no free phase, i.e. any initial disorder drives the polymer into the strong 
disorder phase. In this paper, we will be interested in the low energy excitations above the ground state in dimensions 
d= 1,2,3. 

In disordered systems, there can be states that have an energy very close to the ground state energy but which 
are very different from the ground state in configuration space. For spin glasses, the debate between the droplet 
and replica theories concerns the probabilities and the properties of these states. In the droplet theory ^(|, the low 
temperature physics is described in terms of rare regions with nearly degenerate excitations which appear with a 
probability that decays with a power law of their size. In the replica theory 22], the replica symmetry breaking 
is interpreted as the presence of many pure states in the thermodynamic limit, i.e. the nearly degenerate ground 
states appear with a finite probability for arbitrary large size. More generally, the statistical properties of the nearly 
degenerate excitations (their numbers, their sizes, their geometric properties, the barriers separating them, etc..) 
are interesting in any disordered system, since they govern all properties at very low temperature. In particular, 
a linear behavior in temperature of the specific heat C(T) = bT + 0(T 2 ) seems rather generic for a large class of 
disordered models, including (i) spin glasses where this behavior is measured experimentally [2^ and numerically [24[ 
(ii) disordered elastic systems 25] (iii) one-dimensional spin models where this behavior can be exactly computed via 
the Dyson-Schmidt method 26]. For the last case, the coefficient b of the linear term of the specific heat can be put 
in direct correspondence with the density p{E = 0, 1) of two-level low energy excitations of size I [2?l l28l I29I ] via the 
simple formula b = (ir 2 /6) f dlp(E = 0,1). Other integrals like J dl l k p(E = 0,1) with k = 1,2, ... determine the low 
temperature behavior of other observables. The explicit computation of the density p(E — 0, 1) of excitations as a 
function of their size I has been possible only for one-dimensional models, such as the case of one particle in a random 
potentials [27], and the random field Ising chain via strong disorder renormalization [2^ . For higher dimensional 
disordered systems, the statistics of excitations can only be studied numerically. In particular, there has been a lot of 
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efforts to characterize the distribution and the topology of the low energy excitations in of spin glasses |30| . For the 
directed polymer model in finite dimensions, we are only aware of the work of Tang |31| . where the probability of two 
degenerate non overlapping ground states with binary disorder in 1 + 1 was found to decay as L~ 2 / 3 . Our aim in this 
paper is to measure the statistics pl(E — 0, 1) of low energy excitations of length I, for a directed polymer of length L 
in a Gaussian random potential, in dimensions d = 1,2, 3. We compare our results with the droplet scaling theory in 
finite dimensions We also try to make the connection with the exact results Q on the Cayley tree (d = oo). 

The paper is organized as follows. In section[nJ we briefly recall known results on the directed polymer in a random 
medium, concerning exact results in d = 1 and d = oo (Cayley tree), as well as the spin glass inspired droplet scaling 
theory. After presenting the parameters of our numerical study in section lTTll we briefly present our results on ground 
state energies in section HVI before focusing on low energy excitations. The statistics of boundary and bulk excitations 
are respectively studied in sections IVl and fVll We finally summarize and discuss our results in section IVTT1 



II. SUMMARY OF PREVIOUS WORK ON DIRECTED POLYMERS AT LOW TEMPERATURE 



A. Ground state properties 

The probability distribution of the ground state energy Eq of a directed polymer of length L in dimension 1 + d is 
expected to follow a scaling form 

ft(iihL) ~^( e _*^SS) (i) 

where eo represent the ground state energy density per monomer. The exponent 0^ then governs both the fluctuation 
and the correction to extensivity of the mean value (Note that this is not always the case in disordered systems, see 
e.g. Q]). This result has been proven in d — 1 with the exact value of the exponent [33l 133. 135113^ 

8i = 1/3 (2) 

For the mean-field version on the Cayley tree (d = oo) one has formally 0oo = 02]: the width is of order 0(1), whereas 
the correction to the extensive term Leo m the averaged value En (L) is of order O(lnL) |2|. In finite dimensions 
d = 2, 3, 4, 5, ...the exponent 6d has been numerically measured [3ll l3& l3ii. l4(ij ] . The values of the KPZ exponent \d 
measured in |40j for dimensions d — 2,3 translate into the following values for the directed polymer exponent Od 

2 = 0.244 (3) 

3 = 0.186 (4) 

through the correspondence 9d = Xd/(2 — Xd)- Note that the existence of a finite upper critical dimension has remained 
a very controversial issue between the numerical studies [371 l38l l3& U(| and various theoretical approaches J4l|, |42j, |43J . 
Beyond the exponent 0, the scaling function Vd itself is also of interest : it is exactly known in d = 1 |35l l36| (as well 
as in other geometries 44] ), and has been studied on the Cayley tree, with the conclusion that the distribution is 
not universal but depends on the disorder distribution 45] . Another important property of the ground state in finite 
dimensions d is the probability distribution of its end-point position Rq that follows the scaling form 

Q d (R 0) L)~ 7 ^Qd[? = ^ I ) (5) 



where the exponent Q is directly related to the previous exponent 9d via the simple relation [46| 

Cd = — — (6) 

This corresponds to a superdiffusive behavior £ > 1/2 as soon as > 0, and in particular in one dimension £i = 2/3. 
The probability distribution Qd{Ro, L) has been studied for the case d = 1 in 



B. Exact identities from the statistical tilt symmetry 



In their continuum version, directed polymers in random media belong to a special class of disordered models for 
which exact remarkable identities for thermal fluctuations can be derived 481. The Hamiltonian of these models 
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have a deterministic part which consists in quadratic interactions and a random part whose statistics is translation 
invariant. For the directed polymer in dimension 1 + 1, this so-called statistical tilt symmetry leads in particular to 
the following identity for the fluctuation of the end point 0,010] 

< r\ > - < r L > 2 = TL (7) 

Other identities giving relations between higher cumulants can be similarly derived @, flH l48| . 

The identity Q is rather surprising at first sight, since the fluctuations of the end-point are found to independent 
of the disorder and to be exactly the same as in the absence of disorder! The exact result (0 is particularly interesting 
at very low temperature, since it predicts a linear behavior in T of the second cumulant, and thus puts constraints 
on the statistics of nearly-degenerate excitations. In simpler ID models of one particle in random potentials one can 
explicitly relate [271 l28 j the linear behavior of the position fluctuations to the rare configurations with two nearly 
degenerate minima AE ~ T . For directed polymer in random media, any theory of low energy excitations has to 
reproduce the result Q, and we will now describe its interpretation within the droplet theory @[l|. 



C. Droplet theory for low energy excitations in finite dimensions 

The droplet theory for directed polymers 0, , is very similar to the droplet theory of spin glasses It is a 
scaling theory that can be summarized as follows. The exponent 9 involved in the fluctuations of the ground state 
energy E over the samples (see eq. JTJ) also governs the fluctuations of the energy within one sample as the end 
point varies. As a consequence, if one assumes a scaling distribution analogous to ifT}. the probability to find a 
nearly degenerate ground state AE ~ T of order L behaves as T/L e , so it is rare, but it corresponds to a very large 
fluctuation of the end-point of order Ar ~ with the spatial exponent £ ( see eq. {JjJJ). The contribution of these 
rare nearly degenerate paths to the disorder averag ed fluctuations is thus of order (T/L 6 ) x (L c ) 2 = TL 2C -~ e = TL 
using the scaling relation J§}. 

This naive 'zero-order' argument has to be refined if one is interested into the density pl(E = 0,1) of excitations 
that involve an arbitrary length I for a polymer of length L. For 1 < ! < I, the droplet theory assumes that the 
same scalings apply : the probability of an excitation of length I is of order T/l e and leads to a fluctuation of the 
end-point of order Ar ~ K. However now, the crucial notion of 'independent excitations' [f| has to be introduced 
to obtain a consistent picture. Both for spin glasses ^(| and for directed polymers the idea is that in a given 
sample, droplets with neighbouring sizes tend to have a big overlap. More precisely for the polymer, two excitations 
of lengths and length l<i with ~ I2 will typically merge and then follow the same path to join the ground state, 
as shown by the tree structure of optimal paths to all end-points . As a consequence in [j| , droplets are considered 
to be independent only if there is a factor of order 2 between their sizes, and this gives a factor dlnl = dl/l in all 
integrations over droplets 

dlp md ^(E = 0,l)~^L ( 8 ) 

A more intuitive view of the dl/l factor is to remark that independent excitations stem from a branching process 
along the ground state path, and that the infinitesimal number of branches between I and I + dl is precisely dl/l. 

In summary, the droplet theory predicts a power law distribution of independent excitations, with exponent (1 + 9). 
Note that the absence of any characteristic scale in / means that there exists some infinite correlation length in the 
system in the whole low temperature phase. This is in contrast with simpler models like the random field Ising chain 
|28| or in the spin glass chain in external field p9j | , where the density of excitations was found to decay exponentially 
in /, the correlation length being the Imry-Ma length. 



D. Exact results for low energy excitations on the Cayley tree 

Many exact results for directed polymers on the Cayley tree have been derived by Derrida and Spohn [2] ■ From the 
point of view of excitations, the most important result concerns the distribution of the overlap in the thermodynamic 
limit, which is simply the sum of two delta peaks at q = and q — 1 in the whole low temperature phase : 

n(q) = (l-Y)6(q)+YS(q-l) (9) 

and the distribution of Y over the samples is the same as in the Random Energy Model [49|. In particular, the 
disorder average of eq. © is 0] : 

T „. , / T 



Trfa) = -%)+(!-- )*(?-!) (10) 
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i.e. the overlap is zero with probability T/T c , and one otherwise. This means that for a polymer of length L, the 
important excitations are those of length I ~ L, and that these excitations keep a finite weight in the limit L — > oo. 
To understand the origin of this surprising result, Fisher and Huse Q have computed that the probability to find an 
excitation of length I which branches off at a distance s = L — I <C L from the root behaves as s~ 3 / 2 , i.e. using the 
length notation I = L — s 

p L (E = 0,l)~ {L ^ l)3/2 (11) 



Finally, Tang has studied numerically the overlap distribution Pl (q) for a polymer of finite-size L | 3l| to characterize 
how the two delta peaks develop in i|10[l : the data for < q < 1 follow the scaling behavior (see Figs 4(a) and 4(b) 

of H3) 

^Jqj ~ lT X 'H(q) (12) 
where the scaling function n(q) present the same singularity with exponent 3/2 near q — > and q — > 1 

^ £x(t4^ (l4) 



so that, in the limit L — > oo, ^^(q) only contains two 5 peaks at g = and g = 1 (see eq. I[ll)p). 

III. MODEL AND NUMERICAL DETAILS 

In this paper, we present numerical results for the random-bond version of the directed polymer model on a 1 + d 
hypercubic lattice. The partition function satisfies the following recursion 

2d 

ZL+iif) - ]T e-^ p+e ~>^Z L (r + gj) (15) 

The bond-energies ej / (r'+ gj, r) are random independent variables, drawn with the Gaussian distribution 

*) - ^=e-4 (16) 



For each dimension d, we now give the typical lengths L we have studied, with the corresponding number n s of 
disordered samples : 

(i) For d = 1, L = 50, 100, 200, 300, 400, 600, with respective n s /10 6 = 160, 40, 10, 4.7, 2.5, 4.5. 

(ii) For d = 2, L = 20, 40, 60, 80, 100, with respective n s /10 6 = 30, 4.4, 5.8, 2.2, 1. 

(iii) For d = 3, L = 12, 18, 24, 30, 36, 42, with respective n s /10 6 = 12.5, 2.8, 0.9, 1.6, 0.76, 0.4. 

We first briefly describe our results on the ground state energy statistics, before we turn to the measure of the 
density of low energy excitations on which we focus our attention in this paper. 

IV. DISTRIBUTION OF THE GROUND STATE ENERGY 
A. Scaling distribution 

For the sizes L we have considered, the distribution of the ground state energy Eq follows the scaling form 

as shown on Figs Qfa) , H^b) , Ufa) for d — 1,2,3 respectively. We have checked that the function Fi(u) of Fig. QJa) 
agrees with the numerical tabulation of the exact result given on the web site [5(|. The three functions Fi,F%, F3 are 
shown together for comparison on Fig. |2ft>) . 
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FIG. 1: (a) Rescaled distribution Fi(u) of the d = 1 ground state energy (see equation fT7\ ). for L = 50, 100, 200, 300, 400, 600. 
(b) Rescaled distribution F2(u) of the d = 2 ground state energy (see equation 11711 h for L — 10, 20, 40, 80, 120. 




FIG. 2: (a) Rescaled distribution Js(m) of the d — 3 ground state energy (see equation 1171 '). for L — 6,12,18,24,36. (b) 
Comparison of Ln Fd(u) for d = 1, 2,3 (see equation CZJ). 



B. Behavior of the width AE (L) and average E$ V (L) 



The exponent Od of eq.QJ is expected to govern both the width AEq(L) and the correction to extensivity of the 
average 

AEo(L) ~ L 0d (18) 

^e (d)+i e - 1 e 1 (d) + ... (19) 

Our measures of the exponent 9 a from the width AEq(L) yield 

9 X - 0.33 (20) 

6» 2 - 0.24 (21) 

6» 3 - 0.18 (22) 

are in agreement with the exact value 9% — 1/3 and with the previous numerical measures quoted in Eq. (J1J for 
d=2,3. 
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The fits of the average yields 



d = 


1 : 


E8"(L) 
L 


0.95- 


-0.84 L 01 - 1 


- 1.19 L 


d = 


2 : 


L 


- -1.53- 


- 1.48 L 02 - 1 


- 0.94 L 


d = 


3 : 


E8"(L) 
L 


1.81- 


- 2.05 L 63 ' 1 


- 1.43 L 



(23) 
(24) 
(25) 



V. STATISTICS OF BOUNDARY EXCITATIONS 



Measure of independent excitations 




FIG. 3: Notion of independent boundary excitations (in d = 1) : for each sample, we generate the ground state (bold line) of 
energy Eq, and consider all end-points A% where the energy E(Ai) of the best path ending at A% satisfies E(Ai) — Eq < T (light 
lines). These best paths tend to cluster into families. In our counting procedure of independent excitations, two excitations 
are independent if they have no bond in common. For instance on the Figure, end-points Ai, A2, A3 and A4 count for a single 
excitation associated to the branch point B\. Similarly A7 and As count for a single excitation associated to the branch point 



In this Section, we are interested into the density p b ° undar v ^ — q, Z) of boundary excitations defined as 

1- 



p b ™ ndary (E = 0,l) = hm 



^N h L ° undary {E <T,l) 



(26) 



where AT L ° un ary (E <T,l) is the disorder averaged number of independent boundary excitations of energy < E < 
T — > and of length I existing for a directed polymer of length L. This number A/J" ep (E < T,l) is measured as 
follows. For each sample, we consider all end-points r different from the ground state ro : if the energy E m i n (f) of the 
best path ending at f satisfies E m i n (r) — Eq < T, we construct the best path ending at f to measure its length I, i.e. 
the length over which it is different from the ground state. However, as explained already above when summarizing the 
droplet theory JSJ, one is interested into the number of independent excitations. So here we have used the following 
criterion : two excitations are independent if they have no bond in common. In d = 1, since we use polymers with 
(±) steps, this means that two excitations are independent if they join the ground state at different points (see Fig. 
|3J). In d = 2 and d = 3, two excitations are allowed to join the ground state at the same point provided they branch 
along two different directions. 
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In (|26|l . the parameter T is simply a cut-off and we have checked the independence of the density p boundar v [E — o, I) 
with respect to T for T small enough. For instance in d = 2 and d = 3, we have checked that the cut-offs T = 0.1 
and T = 0.05 yield the same density pboundary ( E = Q, Z). 



B. Boundary excitations for d — 1 




FIG. 4: (a) Log-log plot of the density p^ un aTV (E = 0, 1) of boundary excitations of length I in d = 1 for L = 
50, 100, 200, 300, 400, 600. (b) Same data after the rescaling of equation gJJ with 6i = 1/3. 




On Fig. Ufa), the density p b ° undar y — q, Z) is shown in a log-log plot for various sizes L. These curves can be 
rescaled according to 

p boundary( E = 0, Z) = _L_ R boundary f = *\ ^ 

as shown on Fig. Efb). The master curve R boundar v (x) shown on Fig. Eta) has three important properties 

(i) In the region x — ► 0, the scaling function R boundar v ( x ) follows the power law (see the log- log plot on Fig. Efb)) 

R boundaryM K _J_ ^ 
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so that in the regime 1 <C I <C L, the statistics of independent excitations 

1 

JT+e 



p boundary (E = 0J) ^_l_ ^ j <<; j <<; £ ^ 



follows the droplet power law (JSJ. 

(ii) The function R boundar v jg minimum at some finite value x m i n {d = 1) ~ 0.74, and then grows for x m j n < x < 1 

(iii) In the regime x = l/L — ► 1, the scaling function R boundar v ( x ) diverges. To describe the regime of these very 
large excitations, let us first consider the extreme case I = L of an excitation that branches off at the origin. Our 
result for I = L follow the scaling behavior 

boundary^ = Q [ = £ ^ J_ ^ ^ _ Q.67 (30) 

in agreement with the exponent measured by Tang |3lT ] for the probability of two degenerate ground states in the case 
of binary disorder. More generally, we find that p b ° undar y — q ? I — L — y) with finite y also decays as in (|30|) . see 
Fig. Elb). Since the exponent Ai ~ 0.66 is different from the exponent 1 + 0\ = 4/3 appearing in the scaling form 
(|37Jl . the singularity of the scaling function ^boundary ^ 

near x — * 1 is given by 

R boundary^ K L_ with = 1 + &i — Ai ~ 0.66 (31) 

ir— >i (1 — a;) 171 



C. Boundary excitations for d = 2 




FIG. 6: (a) Log-log plot of the density p b ° undar y(E = 0, Z) of boundary excitations of length I in d = 2 for L = 20, 40, 60, 80, 100. 
(b) Same data after the rescaling of equation l|32|l with #2 ~ 0.24. 



In d — 2, we find the same properties for the statistics of boundary excitations with the appropriate exponent 
6 2 ~ 0.24. The density p b ° undar y( E = 0,1) shown in a log-log plot for various sizes L on Fig. HJa) follow the scaling 
form 

^boundary ( jti n 1\ ^ rjboundary ( „ ^ \ /oo\ 

p y (E = 0, Z) = jj+e;tt \ l) ' 

as shown on Fig. EJb). The master curve n boundar v (x) shown on Fig. Ufa) has the same three important properties 
as in the d = 1 case 

(i) In the region x — > 0, the scaling function follows the power law (see the log-log plot on Fig. 0b)) 

R boundary^ ^ _J_ ,33) 

leading to a statistics of independent excitations 

PL mdarV iE = 0,l)~-L- for 1«Z«L (34) 
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that follows the droplet power law (JHJ. 

(ii) The function R boundar y i s minimum at some finite value x m i n (d = 2) ~ 0.64, and then grows for x m i n < x < 1 

(iii) The density of excitations of length I ~ L that branches off at the origin decays with the power law 

p b ° undarv (E = Q,l = L)~ _L W ith A 2 ~0.53 (35) 
as shown on Fig. Efb). So the singularity of the scaling function R boundar y ( x ) near a; — * 1 is given by 



R boundaryt x \ 1 

V ' X->1 (1 - x) c 



with cr 2 = 1 + 2 - A 2 ~ 0.71 (36) 



D. Boundary excitations for d = 3 




In g? = 3, we find again the same properties for the statistics of boundary excitations with the appropriate exponent 

i ~ 
form 



' 3 ~ 0.18. The density p b ° undar V( E = 0,1) shown in a log-log plot for various sizes L on Fig. Ufa) follow the scaling 



p boundary^ E = Q, |) = _L_ Rboundary ^ = }_\ (37) 
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Ln f(L-y) 




FIG. 9: d = 3 :(a) Rescaled distribution iJ i ™ ndor «(a! = Z/L) (see equation Qfl), that is singular for a; -> and i->1. (b) 
Log-log plot of the density p h ° un ary (E = 0, 1 — L — y) of very large boundary excitations I = L — y with finite y = 0, 1, 2, 3. 



as shown on Fig. IBfb). The master curve R boundar y ( x ) shown on Fig. Eta) has the same three important properties 
as in the previous cases 

(i) In the region x — > 0, the power law (see the log-log plot on Fig. Hfb)) 



^boundary 

leads to a statistics of independent excitations 

boundary / -i-i n t\ 

Pl (E = 0,l) 

that follows the droplet power law (J8J. 



(x) oc — 



+ 03 



for 1< I < L 



(ii) The function R boundar y (#) is minimum at some finite value x m i n {d = 3) ~ 0.56, and then grows for . 



(38) 



(39) 



< x < 1 



(iii) The density of excitations of length I ~ L that branches off at the origin decays with the power law 

boundary^ = Q J^ L ^_±_ ^ ^ _ Q.39 (40) 

as shown on Fig. EJb). So the singularity of the scaling function R boundar y ( x ) near x — > 1 is given by 



I? 



boundary 



(x) OC 



1 



x-»l (1 - x)" 73 



with CT3 = 1 + 03 — A3 



0.79 



(41) 



VI. STATISTICS OF BULK EXCITATIONS 



Measure of independent excitations 



We now turn to the density (E — 0, 1) of bulk excitations defined as 



P ^ lk (E = 0,l) = lim 



(42) 



where Af^ J,lk (E < T, I) is now the disorder averaged number of independent bulk excitations of energy < -E < T — > 
and of length / existing in the bulk for a directed polymer of length L. Here the additional normalisation factor (l/L) 
with respect to the analog definition of boundary excitations <|26|) ensures a coherent normalization between the two 



densities p b ° undar v and 



bulk 



^boundary 
nbulk 



(E — 0,1) represents the probability that the end-monomer belongs to an 



excitation of length whereas p represents the probability that a bulk- monomer belongs to an excitation of length 



/. The number M^ dev {E < T,l) is measured as follows. For each sample, we consider all points fo(t) of the ground 
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-5 1 1 ' 1 ' 1 1 

20 40 60 

FIG. 10: Notion of independent bulk excitations (in d = 1) : for each sample, we generate the ground state (bold line) and 
consider all points Ai of the ground state, with partial energy Eo(Ai). When the energy Ei(Ai) of the second best path ending 
at Ai satisfies Ei(Ai) — Eo(Ai) < T, we construct these bulk excitations (light lines) that join again the ground state at some 
branch point Bj. These bulk excitations tend to cluster into families. In our counting procedure of independent excitations, 
two excitations are independent if they have no bond in common. For instance on the Figure, end-points A4, As, Aq and 
A*; count for a single excitation associated to the branch point B$. The other independent excitations are (^4i,_Bi), (^2,-62), 
(A 3 ,B 3 ) and (A S ,B 4 ). 



state with t — L, L — 1, .., 2 as possible end points of bulk excitations. We consider the best paths joining the ground 
state at fo(t) but arriving from different points than r$(t — 1), to see if they have an relative energy with respect to 
the ground state smaller than the cut-off T . If this is the case, we have found a bulk excitation, and we measure its 
length Z, i.e. the length over which it is different from the ground state. Again, as for boundary excitations, we are 
interested into independent excitations, and we use the criterion according to which two excitations are independent 
if they have no bond in common. In d = 1, since we use polymers with (±) steps, this means that two excitations are 
independent if they join the ground state at different points (see Fig. I10f) . In d = 2 and d — 3, two excitations are 
allowed to join the ground state at the same point provided they branch along two different directions. 



Ln p (1) 
-4 




1+9 
Ln (pL ) 
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FIG. 11: (a) Log-log plot of the density p h ^ lk {E = 0, 1) of boundary excitations of length / in . 
(b) Same data after the rescaling of equation 14311 with exponent 9i = 1/3. 



= 1 for L = 50, 100, 200, 300, 400. 



12 





Lnl Ln(l/L) 

FIG. 13: (a) Log- log plot of the density p b ^ lk {E = 0, 1) of boundary excitations of length I in d = 3 for L = 18, 24, 30, 36, 42. 
(b) Same data after the rescaling of equation l|43|l with exponent #3 = 0.18. 



In dimension d = 1, 2, 3, we find that the density of bulk excitations follows the scaling form 

pb ulk {E = Ql) = -±^R»^ L = A) (43) 

See Figs II 1 1 IT21 1 1 31 for d = 1,2,3 respectively. As d increases, the quality of the rescaling gets weaker, because of the 
smaller sizes L that can be studied via transfer matrix. 

In contrast with the scaling function Jl boundar v ( x ) f boundary excitations, the scaling function R bulk (x) decays 
monotonically for Q < x < 1. In the region x — > 0, the scaling function follows the power law 

R bulk (x) cx (44) 

so that in the regime 1<!<L, the statistics of independent excitations 

p%»<(E = 0,l)~j± gI for 1«/«L (45) 

follows the droplet power law (JSJ). 
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VII. SUMMARY AND CONCLUSIONS 



We have studied the low energy excitations of a directed polymer in a 1 + d random medium. For dimensions 
d = 1, 2, 3, we find that the densities of bulk and boundary excitations follow the scaling behavior p b ™ lk > boundar y = 
0,0 = L- l - e *B! mlk > boundar i>(x = l/L). In the limit x = l/L -> 0, both scaling functions R bulk (x) and R boundar y(x) 
behave as Rb^lk, boundary ( x j _ jg-i-e^ leading to the droplet power law p b ^ lk ' boundar y = 0, i) ~ r 1- *" in the regime 
1 <C I <C I. Beyond their common singularity near a; — ► 0, the two scaling functions R bulk > b ° undar v ( x ) are very 
different (this shows the importance of boundary conditions) : whereas R bulk (x) decays monotonically for < x < 1, 
the function R boundar y ( x ) first decays for < X < .?•„„•„, then grows for .?•„„•„ < .r < 1. and finally presents a 
power law singularity R boundar y ( x ) ~ (1 — x)~ r7d near x — > 1. The density of excitations of length Z ~ L decays as 
boundary^ — QJ = L) ^ L~ Xd where Xd = 1 + 8d — (Jd- Our numerical estimates Ai ~ 0.67, A2 — 0.53 and A3 ~ 0.39 
suggest the relation A^ = 28d, although we are not aware of any simple argument to justify it. However, if it holds, 
this would mean that the scaling function R boundar y (#) has singularities with exponents (1 + 64) and (1 — 0d) near 
x — > and x — ► 1 respectively, i.e. these singularities tend to become the same as 9d decreases, i.e. as the dimension 
d increases (see Fig. El)- This trend is reminiscent of the result on the Cayley tree discussed in Eq. (|12I14|) where 
the singularities have the same exponent on both sides, even if the value (3/2) seems specific to the tree structure 
and cannot be obtained as the limit 9d = in our results. 



R(x) 



20 - 




FIG. 14: Comparison of the scaling function R 



boundary 



x) for d — 1 (thin line), d — 2 (dashed line) and d = 3 (thick line). 



Let us now mention how one recovers the identity J7J of the statistical tilt symmetry. A boundary excitation of 
length / is expected to give rise to a fluctuation of the end-point of order Ar ~ ft (we temporarily drop the dimension 
dependence of the exponents), so that at order T, one gets 



< Ar 2 > - T 



dl l 2 <p b ° undarv (E = 0,l) 



The contribution of excitations of length I with x = l/L finite reads 

-1 



= TL 2C -° / dx x 2<: R boundary (x) 

0<x<l Jo 



< Ar 2 > 

whereas the contribution of very large excitations of length I — L — y with finite y reads 



< Ar 2 > 



TL 2C 



— A 



(46) 



(47) 



(48) 



Since A > 9, the leading contribution is the first one. Using the scaling relation ©, this contribution is of order 
L 2C -~ B = L, as it should to recover |J7J, 

Finally, since the directed polymer model plays the role of a ' baby spin glass', and since various numerical studies 
on spin glasses [3(j find both the droplet scaling behavior for small excitations and system-size excitations governed 
by another 'global' exponent, one may wonder whether both types of excitations can be understood within a single 
scaling function R(x) of the volume fraction x — v/V, where the droplet exponent describes the power law in the 
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regime x — > 0, whereas the statistics of system-size excitations depends on the global properties of the scaling function 
R(x). 
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